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NET-CHARACTEBISTIC  METHOD  PCB  THE  BUMEBICAl  SOLUTION  OP  PROBLEMS  OF 
GAS  DYNAMICS 

K.  M.  Nagoaedov 


Various  nuaerical  aethcds  of  solving  aultidiaensional  probleas 
for  equations  in  partial  derivatives  have  urdergone  intensive 
developaent  and  introduction  in  recent  years.  Different  versions  of 
the  net  aethod  proposed  and  studied  in  the  torks  of  K.  I.  Babenko,  V. 
V.  Rusanov,  N.  N.  Yanenko,  S.  K.  Godunov,  etc.,  the 
Dcrodnitsy n-Belotserkovskiy  aethod  of  integral  relationships,  the 
straight  line  aethod  developed  by  G.  P.  Telenin  et  covorkers,  and 
difference  systeas  using  the  characteristic  properties  of  equations 
[1-9]  are  being  used  successfully  for  solving  specific  probleas  of 
gas  dynaaics. 
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Numerical  set  hods  using  characteristics  have  definite  advantages 
over  ordinary  finite-difference  methods:  in  particular,  the 
consideration  of  the  physical  nature  of  the  prcfclea,  the  small 
variation  in  certain  complexes  from  the  unkcown  functions  along  the 
characteristics,  and  the  possibility  of  predicting  the  tine  of  the 
origination  of  discontinuities,  e.g.,  "suspended”  shocks.  These 
advantages  become  especially  clear  when  solving  problems  with  two 
independent  variables.  However,  the  effectiveness  of  the 
characteristic  method  considerably  decreases  with  the  increase  in  the 
dimensionality  of  space,  and  virtually  vanishes  when  there  are  more 
than  three  independent  variables. 

Thus,  it  becomes  desirable  to  construct  difference  systens  which 
retain  the  main  advantages  of  the  characteristic  method,  but  have  the 
convenience  of  the  net  method  which  results  from  the  fixed  net  and 
the  simplicity  of  the  calculation  formulae. 

Obviously,  the  name  "characteristic  method"  should  only  be  used 
for  numerical  systems  in  which  characteristic  curves  or  surfaces  are 
set  up  during  the  calculation  process.  Then  numerical  systems  in 
which  the  characteristic  relationships  (compatibility  conditions)  are 
only  used  to  derive  the  difference  eguations  at  fixed  nodes  can  be 
considered  to  be  a versicn  cf  the  ret  method.  This  is  even  more 
expedient  because  it  is  precisely  for  this  type  of  system  that  the 
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problem  of  stability  arises,  even  when  the  Courant-Freid richs-Levi 
(CFL)  condition  is  satisfied,  i.e. , the  ranee  cf  the  dependence  for 
the  differential  eguations  lies  inside  the  range  of  the  dependence  of 
the  difference  eguations  [1C-12]. 

There  is  great  flexibility  in  the  selection  of  numerical  methods 
based  on  characteristic  relationships  in  aultidinensiona 1 problems. 
Report  [17]  gives  a rather  detailed  survey  cf  the  specific  methods  of 
this  type  which  have  been  developed  for  gas-dynamic  eguations.  But, 
obviously,  it  makes  sense  to  have  a sufficiently  general  and  simple 
method  of  constructing  explicit  difference  systems  for  hyperbolic 
multidimensional  eguations  similar  to  the  Ccurant,  Isacson,  Rees 
method  [13]  or  the  "running  calculation"  method  (see  [15])  for 
systems  with  two  independent  variables. 

The  main  idea  of  the  net-characteristic  method  proposed  in  [|8  ] 
is  that  difference  approximations  of  the  ccmpatibility  conditions  and 
linear  or  guadratic  interpolation  by  the  nodes  of  a preassigned  fixed 
network  are  used  for  each  elementary  cell.  Cn  the  other  hand,  the 
method  can  be  considered  to  be  the  ordinary  method  of  nets  with  a 
specific  rule  for  the  difference  approximation  cf  the  initial 
equations.  This  is  the  viewpoint  from  which  the  net-characteristic 
method  of  constructing  difference  systems  is  considered  below. 
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§1  gives  the  sain  concept  of  the  net-ckaracteristic  Method  and 
soae  forms  of  coapact  notation  of  difference  equations  of  the  initial 

j 

Matrices  using  fundaaental  aatrices.  §2  investigates  the  local 
stability  cf  the  corresponding  difference  systeas.  §3  shows  that  the 
net-characteristic  Method  can  be  used  to  cccstruct  divergent 
difference  systeas.  First-order  systeas  are  the  siaplest,  being 
acnotonic  and  flexible. 


A number  of  probleas  of  gas  dynamics  have  keen  solved  by  the 
net-characteristic  method: 

- supersonic  steady  flew  about  both  ccnvex  and  concave  sharp  and 
bluff  bodies  at  an  angle  of  attach; 

- the  direct  problea  of  calculating  flow  in  a nozzle; 


- three-diaensionai  flow  about  segaented  ccnical  bodies,  etc. 


The  last  probleas  were  solved  by  the  deter a inat ion  method  which,  in 
the  aost  interesting  cases,  aeans  solving  equations  in  partial 
derivatives  with  four  independent  variables.  §4  gives  soae  results  of 
the  calculations  as  an  illustration. 


5 


§1.  Net-Characteristic  Method  of  Constructing  Difference  Systens  for 
Cuasilinear  Hyperbolic  Equations 


First  we  will  consider  a systea  of  first-crder  hyperbolic 
guasilinear  equations  in  partial  derivatives  with  two  independent 


variables  t,  x 


1i«*  *U.-i 


where  the  aatrix  A and  the  vector  colunn  of  the  right  sides  f can  be 
functions  of  t , x and  vcctot  function  u of  tnknown  variables  with  N 
coaponents. 


Let  and  u • (i  * 1,  2,  ...»  N)  be  prcper  real  nuabers  and  the 
independent  latent  vectors  cf  transposed  aatrix  K*  corresponding  to 
then.  Multiplying  (1.1)  in  scalar  fora  by  u-(  and  considering  the 
equation  , we  will  reduce  the  systea  to  the 

characteristic  (normal)  fort 


to,  vt*. - u). I *>)  a.t) 


Suppose  that  we  knew  solution  (1.1)  in  layer  t = t0  •*  const.  He  will 
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find  the  value  of  0H  * U ( H)  at  point  H of  a close  layer  t = t0  ♦ r. 
Approximating  (1.2)  by  difference  equations  at  , ve  will 

have 

Here  the  subscript  i designates  the  values  cf  the  functions  at  the 
points  of  contact  of  the  characteristics  drawn  fron  point  H with  line 
t * tD.  He  will  point  out  that  it  is  necessary  to  coapute  «j  fro*  the 
aean  parameters  of  points  H and  i in  order  to  oktain  the  second-order 
approximation  for  quasilinear  equations  at  v * l/2  in  (1.3). 

Based  on  the  linear  independence  of  vectors  0^,  systea  (1.3)  has 
a single  solution  relative  to 

11m * p (t,  x,  u(t»,  x),  uH,g)*0r«-2m*» 

He  will  consider  the  fixed  network  ), 

x«  mM'n*  o,44,  *2,..  ) and  ve  will  designate  the  values  of  the  functions 
at  the  nodes  of  the  network  by  11*  . If  the  relationship  of  the 
paraaeters  with  subscript  i to  the  known  values  of  ti*  in  layers  t = 
tn  * nr  is  indicated,  we  can  find  a difference  systea  with  a fixed 
network,  or  the  net  aethod.  When  using  three  points  for  linear  or 
quadratic  interpolation,  we  will  have,  respectively  (*  = r/h) : 
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r>  n 


tC-o>.  gu%o/ 

Relationships  (1.3)  and  (1.4)  completely  detine  the  difference 
system,  but  it  is  expedient  tc  consider  its  compact  form.  We  will 
introduce  the  following  designations 


*.u.v«<t 0J£\ 

ti  - y*  (u z..  - ui_,  - £uZ) 


Obviously,  matrix  Q,  whose  lines  are  the  linearity  independent 
vectors  wj,  is  not  confluent  and  has  the  inverse  matrix  Q"1.  Writing 
(1.3)  in  matrix  fora  and  using  the  result  ot  (1.4), 

Q 8u;-«ACc»u*  s'a'Q^  u 

where  y = 1 (2)  in  the  case  of  linear  (quadratic)  interpolation, 

A . A*  are  diagonal  matrices  with  elements  { u,.  ./•••*}  and 

Iju.l*. ...  lyL.i* J , we  can  find 
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tC--  lC- 5 QHAQ*m  a 4 s'S2*AQ.':m  at  u.. 

♦ t{W^-tl-v)C‘-  VS’fl’AOU}* 

The  compatibility  conditions  (1.2)  in  a fixed  system  of  nodes  using 
interpolation  formula  (1.9)  were  used  to  derive  difference  systems 
(1.5).  But  the  difference  formulae  can  be  written  immediately  if  we 
know  the  fundamental  matrix  Q for  A,  i.e.  , the  matrix  which  reduces  A 
to  diagonal  form.  We  will  pcint  out  the  two  most  important  systems  in 
(1-5):  first-order  system  I with  linear  interpolation  (y  = 1,  v = 0, 
then  A-Q'xaQ  and  f can  be  computed  at  pcint  t = nr,  x = mh 
without  disturbing  the  crder  of  the  approximaticn) , and  system  II 
with  quadratic  interpolation  (y  = 2,  \)  - »/2)  „ 

We  will  now  consider  the  multidimensional  case 
U«+  4,uXfw..-f  u.fc5 

we  will  limit  ourselves  to  the  terms  described  in  (1.6)  when  deriving 
difference  formulae  of  the  type  in  (1.5),  since  the  extension  of  the 
method  to  a large  number  of  variables  is  dene  formally.  We  will 
assume  that  matrices  X(f  A2,  ...  have  only  proper  real  values.  Let  ju* 
and  u*  (i  * 1,  2,  ...,  N)  be  proper  numbers  and  the  independent 
latent  vectors  of  matrix  A»  (a  = 1,  2)  corresponding  to  them.  Then 
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equations  (1.6)  can  be  reduced  to  cne  of  the  following  characteristic 
f oris 

L>!  U,  - U>i  ( f - A,Ui,t 

u*t,  *) 

H u*.i  * (* ” “*■* J 

r “itj  * Ul  » ^ ti  t,  ; j * i.  2.  ■ • •>  *1) 

Here  U»,  , uti  designate  the  derivatives  along  the  two-dimensional 
character! sties  on  the  surfaces  x2  = const  and  x,  = const, 
respectively,  or,  in  other  words,  along  the  lines  of  intersection  of 
the  characteristic  and  coordinate  surfaces.  The  directions  of  these 
lines  are  deterained  by  the  vectors  and  {i.o./u,1} 

He  will  now  consider  3N  equations  (1.6)  — (1.7)  -jointly. 

Obviously,  the  nuaber  of  independent  equations  among  them  is  equal  to 
N.  This  fact  can  be  used  to  construct  converient  explicit  difference 
systems  for  multidimensional  equations,  in  (articular,  those  not 
requiring  the  approximation  of  the  partial  derivatives  for  Xi  and  x2. 
Replacing  equations  (1.6)- (1.7)  by  the  finite-difference 
relationships  along  the  corresponding  lines  at  t- t,*U-v)Z  and 
multiplying  (1.7)  by  (2*  is  a na  trix  whose  lines  are  vectors 

( oo*  (<*•*,*  ) , we  will  have: 


1 
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A4  UxJh+  \'[{-  A,  Uxt-A£UZe)*- 


Un-  S'/tQ.tU) 


* (1-  VJC'C - **Ux4l„  - A|tlx£  ]/• 


Subtracting  the  first  equation  fron  the  sub  of  the  last  equations  and 
using  interpolation  foriulae  (1.4)  to  find  £,u  and  Q,  U;  we  will 
obtain  the  unknown  foraula 


tC  - U«-%  «!aC,  a»  u ♦ 5?  2<  A*  Q At  11  - 

- 6,  Ci’A,  Q,  Ac  U + Sf  q?a[  C At  u + 

♦ t {vC  *a-V)C‘-  VfS,  A**t  fj  + 

+ VC*T41(A£Ux1)x<  ♦ Aj(^UOJt4  0ra-£V)tl-*E*J  (i.8) 

Here  1,2) , I,*  "•K  5 ij-£h*.then  deriving  (1.8),  we  consider  that  it 

suffices  to  take  the  coefticieat  at  Vr*  with  accuracy  down  to  the 
terns  cn  the  order  of  0 (r) . Furthetnorc,  this  term  sust  only  be 
considered  for  seccnd-oider  system  II  J,  V =*»/*)#  whereupon  it 

is  necessary  to  use  one  more  point  in  layer  t = nr  to  calculate  the 
mixed  derivatives. 


f 
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It  should  also  be  pointed  out  that  for  the  first-order  system 
(•\>-  0)  - system  I - analogously  tc  (1-5)  and  (1.8),  we  can 
immediately  write  the  formulae  for  any  number  of  three-dimensional 
coordinates,  as  long  as  natrix  A*,  is  reduced  to  the  diagonal  form 

(aC.J.X,  ...) 

§2.  Stability  of  Net-Characteristic  Method 

He  will  consider  the  differential  eguatiens  (1.1)  and  the 
difference  equations  approximating  them  (1.5)  at  A = const,  f * 0 and 
periodic  boundary  conditions.  Then  the  stability  and,  therefore,  the 
convergence  of  the  difference  systems,  which  follows  from  it 
according  to  the  to  the  equivalence  theorem,  can  be  obtained  by 
analyzing  solution  (1.5)  in  the  form  [14] 

U-m  c G?U°  ( f • kK,  0,  .. .) 

G * 5 (C,  «f)*E-5rQJ.\*Q  (i-coStp)-;6S2VkQ 


Neumann's  necessary  condition  of  stability 

mo.x  | A,|«  max  {ri-fi’l/ijVl-eosipVf \ 


COC  = 0951 


PAGE  12 


where  X ■ is  the  proper  value  of  the  transitional  matrix  G,  will  be 
satisfied  at  ffl/t.ui  (i  = 1,  2,  ...,  N)  , i.e.,  when  the  CFL 
condition  is  satisfied.  This  sane  condition  is  also  sufficient,  since 
liaitation  Gn,  unifora  with  respect  to  t,  fellows  from  the 
relationship 


Thus,  systems  I and  II  of  the  net-characteristic  method  are 
stable  in  the  case  of  two  independent  variables.  We  will  point  out 
that  for  the  first-order  systea  - system  I - the  stability  of  the 
Cauchy  difference  problem  can  be  proven  with  less  stringent 
liaitations.  The  proof  of  the  convergence  of  a system  similar  to 
(1.5)  is  given  in  [13]  for  quasilinear  system  (1.1).  At  j - 1, 
relationship  (1.5)  can  be  written  as 


tU  ■ T,  ♦ T,  tu.  +TtuZ..  ( ) 


Since  T,  (i*0,ti)  ^e  matrices  nonnegative ly  determined  at 


£ HA«*i 


, (2.1)  is  a positive  Friedrichs  system.  Book  [15] 


considers  certain  similar  systems;  in  particular,  the  running 
calculation  system,  which  is  similar  to  systfm  (1.6)  or  systea  [13]. 


A 
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He  will  now  consider  the  case  of  three  independent  variables 
(1.6).  At  V * 0 ( y = 1),  for  (1-8)  the  traositional  matrix 

S«  6 (t,*.  ,f,)  ( H*,  ..,+**.2) 

has  the  forn: 

G = E - 6i(i-cob\f1)A,  - c 6^  bin.  'PiA,~ 

- ff*  (l-costf,)  wa«fvA*. 

A.  = A1.  Q*  (<*=1.2) 

He  will  find  the  Neumann  stability  condition  with  the  assumption  of 
the  symmetry  of  matrix  4„.  Let  z * x ♦ iy  be  the  latent  vector  G. 
Considering  equation  (Gx.I)- A |zl*  and  the  symmetry  of  matrix  At,  we 
can  find 


p _ C4«  x,  x ) ♦ A.  V V) 


| 


4 fls,  liJVl  * Qec  (<*»!.*) 


1 
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It  is  evident  fro*  inequality  (2.2)  that  the  Neuaann  stability 
condition  will  be  valid  at  s,  a , ♦ \ CL,  < »*»a(s,iju‘i+e,ifj*|)fc. i,i.e.  , 

vhen  the  CPL  condition  is  satisfied.  The  syseetry  of  A was  used  in 
deriving  this  equation.  But  the  considers ticn  of  individual  examples 
of  asyssetric.al  aatrices  wakes  it  possible  tc  hypothesize  that  the 
Neumann  condition  is  satisfied  for  system  I (1.6)  when  the  CFL 
condition  is  satisfied,  and  in  general. 

In  the  case  of  one  scalar  equation  (1.6)  (N  = 1),  stability 
condition  (2.2)  assumes  the  form 

)A|*«  f 1 /I*  ( i - co*  ) - 6k  /!«  ( 1-  cos  <fk  )j'-t 

+ Cff4,(jsi«v  *,1  «•  ff,  A,  ismifcl]  4 * i (2.3) 

This  relationship,  which  is  a necessary  and  sufficient  condition  for 
stability  in  the  case  of  corstant  coefficients,  is  equivalent  to 
condition  (2.2),  at  least,  for  symmetrical  aatrices.  Therefore,  the 


1 


I 
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■ain  properties  of  the  difference  systems  of  the  method  proposed 
above  can  be  revealed  in  the  simplest  example 

A*ttV  ( A**0,  ot*i.a;  (2.4) 

Adding  a new  point  (■  - 1,  / - 1),  selected  with  consideration  of  the 
direction  of  the  characteristic,  tc  the  five  symmetrical  points  in 
layer  t -t*  , the  difference  systems  of  net-characteristic  method 
(1.8)  can  be  written  as  fellows  for  equaticr.  (2.4): 


At  V »0,  t * 1 we  have  systea  1.  It  is  easy  to  check  the 
stability  of  this  system  when  the  CFL  condition  is  satisfied  by  the 
subscript  of  the  difference  system,  since  t)«  system  is  positive: 

I - I i-€Li -o«l ♦ a, ♦a**  cc  u !♦  ct 

*1 , c-  ma*(d 

he  will  point  out  that  this  system  is  mcnotcnic,  i.e.,  the  aonotonic 
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profiles  are  also  transformed  into  monctonic  profiles  at  fs0 

Using  the  Fourier  method,  i.e.,  relationships  like  (2.3),  we  can 
study  the  difference  systems  obtained  from  (2.5).  We  will  point  out 
seme  of  them,  omitting  the  calculations. 

During  quadratic  interpolation  (y  = 2) , the  system  turns  out  to 
be  stable  when  the  CFL  condition  is  only  satisfied  for  V = ‘/a 
(system  II).  The  system  is  stable  for  the  rest  of  the  values  of  p’ 

when  a*  = 0(r),  i.e.,  when  the  tine  interval  is  considerably  smaller 

< 

than  the  spatial  coordinate  interval.  Then  i + a*+  a\  , and  the 

growth  of  the  initial  errer  i0  during  the  calculation  is 
characterized  by  the  value  £.(!•» a* 4ai  )** 

The  use  of  additional  nodes  on  plane  t * nr  lowers  stability. 

For  example,  thy  use  of  linear  interpolation  with  the  satisfaction  of 
the  CFL  condition  when  using  yet  another  point  (a  - 2,  i ) makes  the 
system  generally  unstable.  Systems  which  arc  only  stable  at  a'  - 0 (r) 

I 

are  obtained  during  high-order  interpolat  icc  using  new  points. 

Systems  can  exist  which  are  stable  when  the  CFL  condition  is 
satisfied  in  one  direction,  e.g.,  a,  < 1,  and  a*  * 0(v)  is  required 
for  the  other  variable.  Specifically,  this  is  the  system  previously 
proposed  by  the  author  of  [8].  During  the  practical  calculation  using 
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this  systea,  the  interval  ot  variable  x2  was  selected  on  the  basis  of 
the  nature  of  the  solution  cf  the  problem  sc  that  a*  - 10"*.  This 
results  in  a 12o/o  increase  in  the  initial  error  f0  after  1000 
intervals  on  t,  which  is  completely  permissible.  This  example 
indicates  that  these  systeas  can  be  used  when  a*  <<  1. 

The  nature  of  the  suppression  (or  propagation)  of  the  errors 
arising  during  the  calculation  is  an  iaportant  property  of  any 
difference  systea,  as  well  as  the  behavior  cf  the  net  functions  in 
the  region  of  an  irregular  solution.  Rela t ienships  (2.5)  and  the 
practical  calculations  show  that  systea  I acnotcoically  decreases  the 
error  (this  is  logical,  since  this  systea  is  positive),  while  systea 
II  also  decreases  the  aaplitude  cf  the  error,  but  with  fluctuations 
cf  different  signs.  A similar  situation  occurs  for  aore  complex 
equations.  Here  it  should  be  pointed  out  that  this  behavior  of  the 
soluticn  cannot  be  considered  to  be  instability,  since  se  can 
strictly  prove  stability  in  (2.5). 

Based  on  the  above  discussion  and  the  analysis  of  the  acoustic 
equations  made  in  [ IB],  we  can  draw  the  following  conclusions. 

1.  From  the  standpoint  of  stability,  the  best  systea  is  the 
first-order  systea  with  linear  interpolaticn  - systea  1.  This  systea 
is  aonotonic  and  contains  the  dissipative  terns  necessary  for  finding 
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1 


solutions  which  are  not  very  regular. 

2.  Systee  II  can  cause  unfavorable  fluctuations  for  irregular 
solutions,  but  nevertheless  Bakes  it  possible  tc  use  larger 
intervals.  Therefore,  it  is  advisable  to  use  this  systen  for 
nuaerically  determining  rather  regular  solutions. 


I 


3.  It  is  expedient  tc  write  a single  ccmputer  progran  for  the 
net-characteristic  nethcd,  and  select  the  order  of  interpolation  on 
the  basis  of  the  nature  of  the  pr  tlem. 

V Using  quadratic  (or  higher-order)  interpolation  at  V*  i/i 
Bakes  the  difference  systen  unstable  at  a;  = const  and  stable  only  at 
a?  * 0 (r) . Checking  for  many  exaaples  of  hyperbolic  equations  showed 
that  the  first-order  systews  for  t and  higher-order  systems  for  xt , 
x,  are  only  stable  at  a*  = 0 (r) . 

: i 

I i 

§3.  Divergent  Systems  of  the  Net-characteristic  Method 

The  concept  of  the  net-characteristic  method  can  also  be  used  to 
construct  divergent  difference  systems.  First  we  will  consider  a 
systen  of  hyperbolic  equations  in  divergent  form  with  two  independent 
variables 

ii J 
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Ut  + Fx-{  F*F(t,  X.U.),  {••((t.x.u)  (3.1) 

Suppose  a Matrix  with  fixed  components  n,  assigned  in  space  U,  has 
the  peeper  values  p(i  = 1,  2,  ...»  N)  with  the  sane  sign  as  the 
proper  values  of  matrix  Fu#  i.e.,  they  have  the  same  diagonal  matrix 
L» {«*«"/•», Sign/%v..,  kign/*.}  • Hultiplying  (3.1)  ty  linearily  independent 
latent  vectors  w(  of  matrix  H in  scalar  fort,  we  will  have 

<*>» wt  Fx - co4  f ((.«  W)  (3.2) 

Based  on  the  independence  of  u^,  this  systex  is  equivalent  to  the 
initial  system  (3.1).  Approximating  (3.2)  with  consideration  of  the 
characteristic  directions  of  matrix  H by  tbe  one-sided  differences 
relative  to  x,  we  will  have 


uT-  ui 


U>: 


F - C r 


c"  - F 

» m in 


at 

at  >o 


: oJ 


l* 


CO, 


tC-  u" 


. ur"  ~ u-f.  UJ.  r _ . . 

1 C + h.  [ A">F~  S.gn/U^Ap,  F]=u>t{ 


L_ 


Hultiplying  the  difference  eguatiens  on  the  left  written  in  matrix 
form  by  Q-*  and  using  the  previously  introduced  designations,  we  will 
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obtain  the  first-order  divergent  systea 

C*u:-6a.F  + ffQ-‘LCA!.F^C 

He  can  consider  the  difference  system  for  version  U in  order  to  study 
the  local  stability  of  (3.3),  and  the  systea  will  be  stable,  e.g.,  at 
H ■ I * Fu.  Daring  practical  calculations,  for  each  point  (n  ♦ 1,  a) 
a#  can  use  or  retain  it  for  other  ■ when  the  proper  values  of 

Fu  do  not  change  sign. 

For  systeas  of  hyperbolic  eguations  with  a large  nunber  of 
variables 

+ <F«)x, * i (3,4) 

the  first-order  divergent  systea  can  be  writtec  analogously  to  (3.3) 
and  (1.8): 
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U*  - - e4  Ft  ♦ 6.  S?;‘  Lj  c, At  Fi- 

“ ®<e  A«  p«  ♦ 6,  C;‘L4Q  Aj  R,  4- 


+ c*^  -t-  0(T+  hUj'f  h*) 


We  will  point  out  that  the  simplest  fota  of  divergent  systems  is  used 
when  solving  problems  in  the  regions  in  which  the  proper  values  of 

the  matrices  (F,)u  . (F»)u do  not  change  sign.  Then  it  suffices  to 

calculate  the  matrices  in  (3.3)  and  (3.5)  at  one  point. 


We  will  compare  the  divergent  system  of  the  net-characteristic 
ithod  with  the  Laics  system.  For  equation  (j.  1) , the  Lais  system  is: 


A i ^ ® ^ — * 1 ^ I a r 1 i ^ - r~  l . . _ alt 

11  „ * r O Awe  = lim 


This  system  is  stable  when  the  CPL  conditicc  is  satisfied,  but  unlike 
(3.3)  , it  is  not  flexible,  i.e.,  fcr  small  «,  which  is  unavoidable 
with  nonlinear  equations  due  to  the  selecticn  cf  the  minimum  interval 
using  the  CFL  condition,  it  approximates  a parabolic  equation  (see 


1 
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[15]).  The  absence  of  • in  the  seccnd  diffeience  in  (3.6)  also  causes 
the  branches  to  be  "spread"  into  a larger  number  of  steps  than  in 
systea  (3.  3)  . 

In  conclusion*  we  will  give  the  second-order  divergent  system 
for  equation  (3.4)  without  deriving  it: 

*1  I A • 

*?»]■•■  ora**  *;.►  i,; ) 

U.”> 

^»  = (F»)u  , fla  = (Fa)u 

CP  = F,)*^  - A,  ( V ( + 

The  values  of  Al#  A2  and cp  are  calculated  at  point  ( n,  a,  l ) with 
precision  down  to  the  terns  on  the  order  of  0 ( h t *h2). 


r 
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Gas  Flews 

The  above  method  can  be  used  to  study  the  movement  of  an 
inviscid  non-heat- conductive  gas  ic  the  case  of  purely  supersonic 
steady  flows,  as  well  as  in  the  general  case  of  unstationary 
gas-dynamic  equations.  The  corresponding  difference  equations  of  the 
net-characteristic  (NC)  [CX]  method  are  partially  given  in  [18], 

[19],  while  the  method  of  calculating  the  boundary  points  (on  the 
surface  of  the  body  and'  cm  the  shock  wave)  is  described  in  [8],  [15]. 

Only  some  of  the  most  typical  examples  cf  the  numerical  calculations 
are  given  below. 

The  main  properties  of  the  NC  method  were  studied  in  the  problem 
of  supersonic  steady  flew  cn  the  side  surface  of  blunt  cones.  Many 
results  were  gathered  for  this  problem,  and  different  numerical 
methods  were  obtained  [20].  The  calculations  showed  that  the 
parameters  on  the  bedy  and  the  shock  wave  are  computed  with  an  error 
of  1-2o/o  both  using  system  I,  and  using  system  II  of  the  NC  method. 

But  here  system  I requires  approximately  four  times  as  many  nodes  as 
system  II  for  obtaining  the  same  degree  of  accuracy.  But  the 
seccnd-order  system  results  in  fluctuations  ("ripples")  around  a 
certain  averaged  curve  which  are  especially  noticeable  when 
calculating  the  profiles  cf  gas-dytanic  functicns  in  regions  with  a 

narked  change  in  parameters  and  with  discontinuities  in  the 

I 

L L 
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derivatives.  The  pattern  shewn  in  Fig.  1 is  the  most  characteristic. 
This  figure  shows  the  pressure  profiles  P cn  different  meridional 
planes  if  = const  at  a distance  of  five  radii  of  blunting  for  a 
spherically  blunted  cone  with  a half-angle  cf  aperture  of  10° 
streaalined  by  a supersonic  gas  flew  with  an  adiabatic  index  of  t 5 
1.4  at  a Hach  number  of  It  = 6 and  an  angle  cf  attack  of  a = 5°.  The 
solid  lines  show  the  results  of  the  calculation  using  the  direct 
systea  of  the  characteristic  method  [21],  while  the  crosses  and 
circles  indicate  the  data  calculated  using  systems  I and  II  of  the  NC 
method.  It  is  evident  that  while  method  [2  1]  clearly  isolates  the 
discontinuities  (they  originate  due  to  the  interruption  of  curvature 
at  the  point  of  contact  of  the  cone  with  the  sphere),  the  NC  method 
smoothes  these  discontinuities , whereupon  system  I is  regular,  and 
system  II  has  marked  fluctuations. 


■ 
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A similar  pattern  is  also  observed  when  performing  a "through" 
calculation  of  the  zones  of  discontinuities  in  the  functions.  In 
ctder  to  test  the  possibility  of  applying  tie  KC  method  to  the 
through  calculation  of  disccnt inuous  solutions,  various  problems  were 
considered.  Figure  2 shows  the  flow  pattern  about  a two-degree  cone 
at  II  = 3 and  j = 1.4,  as  well  as  the  pressure  distribution  P on  the 
body.  The  internal  connections  of  the  shocks  were  found  by  the 
through  calculation  using  system  I.  The  solid  lines  show  the  results 
frcn  [22].  Figure  3 shows  the  flow  pattern  about  a plane  body  and  the 
pressure  distribution  on  the  body  at  H = 5.  Here  the  leading  shock 


wave  1 was  isolated  by  calculation,  while  tie  "suspended"  shock  wave 
2 was  obtained  by  the  through  calculation.  The  upper  curve  2 depicts 
the  precise  position  of  the  shock. 
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The  following  should  be  pointed  out  regarding  the  use  of  the  NC 
aethod  for  a through  calculation.  The  calculations  shoved  that 
discontinuities  which  ace  not  very  narked  can  te  obtained  using 
systea  I of  the  NC  aethod.  This  system  is  net  divergent,  but.  in  the 
difference  systea,  the  pressure  and  angles  cf  slope  of  the  velocity 
vector  for  steady  flows  [18]  are  deternined  independent  of  the  energy 
equation.  Here  the  coefficients  in  these  equations  depend  mainly  on 
the  local  Mach  nuaber,  and  vary  little  in  many  problems.  The  energy 
equation  is  used  in  a such  a fora  that  the  dissipative  term  in  the 
presence  of  a shock  wave  results  ir  the  Hugcnict  condition  for  the 
relationship  between  pressure  and  enthalpy.  If  it  is  necessary  to 
cross  marked  discontinuities,  e.g.,  a leading  shock  wave,  by  a 
through  calculation,  it  is  necessary  to  use  the  divergent  systea  of 
the  NC  aethod  given  in  §3  tc  increase  accuracy. 

As  we  know,  aany  problems  of  gas  dynaaics  require  the 
siaultaneous  calculation  of  the  regions  of  subsonic  and  supersonic 
flow.  In  order  to  make  it  possible  to  apply  the  NC  aethod  to  these 
probleas,  we  can  introduce  a new  independent  variable  t such  that  the 
equation  becomes  hyperbolic  everywhere  and  the  solution  at  t -*  - gives 
us  the  solution  of  the  initial  equations.  This  is  done  siaply  by 
retaining,  e.g.,  physical  tiae  t,  tor  gas-dynaaic  equations.  This  is 

i 
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precisely  how  problems  of  three-dimensional  flew  about  segmented 
conical  bodies  and  the  direct  problem  of  flew  ir  a nozzle  were  solved 
by  the  net-characteristics  method. 

Figure  4 gives  the  flow  pattern  for  a segmented  conical  body 
shaped  like  the  Apollo  launch  capsule,  which  is  streamlined  at  an 
angle  of  attack  of  a = 0.  The  shock  waves  and  scund  lines  (a)  are 
shewn,  along  with  the  pressure  distribution  on  the  body  (b)  for  a 
perfect  gas  at  flach  numbers  of  the  incoming  flew  of  fl  = 2,  6 and  20 
and  adiabatic  indices  y - 1.4,  as  well  as  with  consideration  of  the 
equilibrium  physico-chemical  processes  in  the  air  at  velocities  of 
the  incoming  flow  of  V,  » 5 km/s  (FI  = 16.1)  and  V-  * 7.5  km/s  (H  = 
23.5).  Similar  results  are  shown  in  Fig.  5 at  B * 6 and  y = 1.4  for 
three-dimensional  streailining  at  angles  of  attack  of  a = 10°,  20° 
and  25®. 


J 
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Pl9.  4. 


III! 
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Figure  6 shows  another  problea  - the  calculation  of  a circular 
nozzle  with  an  assigned  shape.  It  gives  the  shape  of  the  nozzle,  the 
sound  line,  and  the  pressure  distribution  and  longitudinal  velocity 
cn  the  upper  and  lower  walls. 


Fig.  6. 
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Conclusion 

The  net-characteristic  nethod  proposed  for  solving 
aultidiaensional  hyperbolic  equations  lakes  it  possible  to  calculate 
different  gas-dynaiics  problems.  This  report  focused  nainly  on  the 
properties  and  peculiarities  of  the  difference  systea,  and  the 
results  given  are  illustrative,  but  they  nevertheless  lake  it 
possible  to  deteraine  the  class  of  probleas  which  the  system  can  be 
used  to  solve.  The  analysis  and  nunerous  experiiental  calculations 
show  that  the  method  is  stable  and  ncnctonic,  and  that  it  makes  it 
possible  to  perfora  "through"  calculations  which  do  not  generate 
peculiarit ies. 

This  aethod  was  developed  at  the  Bcsccw  Physi co-Technical 
Institute  as  part  of  the  numerical  aethod  research  plan  conducted 
under  the  guidance  of  Professor  0.  fl.  Belot rerkovskiy.  The  aain 
concepts  of  the  aethod  are  published  in  the  report  by  K.  R.  flagoaedov 
and  A.  S.  Kholodov  [18].  A.  S.  Kholodov,  V.  I.  Kosarev  and  V.  V. 
Pirogov  conducted  aany  of  calculations  using  the  net-characteristic 


aethod 


I 
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